function [ A_ineq,b_ineq] = constraint_mat_theta_2(combos_all_unique,J)

% monotonicity

[ A_monot ] = A_monot_fun(combos_all_unique,J);

A_ineq = sparse(A_monot);

b_ineq = [-ones(size(A_ineq,1),1)*0.0];

end

